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ABSTRACT 



This report illustrates the numerical procedure used by the Aero- 
Astrodynamics Laboratory in performing roll stability analyses for Saturn type 
space vehicles. Vehicle control in the roll channel is provided by a position 
gyro and a rate gyro. The dynamics of the rate gyro are included as a separate 
equation in the equations of motion. The modes of oscillation in this system 
are (1) torsion, (2) roll, (3) roll sloshing, (4) rate gyro roll, and (5) swivel 
engine. Two methods are available for solving the system's eigenvalues: 
the matrix iteration method and the characteristic equation approach. In 
both cases the system is solved for the normalized eigenvectors, which pro- 
vide additional information about the system. 


The frequency response of the system, at the control thrust vector 
point, can be found for a sinusoidal forcing function at various sensor lo- 
cations. 


The equations of motion are derived in the appendix. 
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TECHNICAL MEMORANDUM X-53368 


NUMERICAL PROCEDURES FOR ROLL STABILITY STUDIES 

SUMMARY 


This report illustrates the numerical procedures used by the Aero- 
Astrodynamics Laboratory in performing roll stability analyses for Saturn type 
space vehicles. Vehicle control in the roll channel is provided by a position 
gyro and a rate gyro. The dynamics of the rate gyro are included as a separate 
equation in the equations of motion. The modes of oscillation in the system are 
(1) torsion, (2) roll, (3) roll sloshing, (4) rate gyro roll, and (5) swivel 
engine. Two methods are available for solving the system’s eigenvalues: the 
matrix iteration method and the characteristic equation approach. In both 
cases the system is solved for the normalized eigenvectors, which provide 
additional information about the system. 

The frequency response of the system, at the control thrust vector point, 
can be found for a sinusoidal forcing function at various sensor locations. 

The equations of motion are derived in the appendix. 


INTRODUCTION 


The stability analysis of an elastic space vehicle requires the eigenvalues 
for the solution of a set of linear, homogeneous differential equations. These 
differential equations describe the dynamics of the space vehicle and the 
characteristics of the control system. 

The solution is readily obtainable for a simple system, but it is difficult 
for systems with a more complex control system and large numbers of oscilla- 
ting modes because of the high order of the characteristic equation and the large 
magnitude of the coefficients. 

To adequately describe the vehicle characteristics, the following modes 
of oscillation are included in the system: 



torsion 

<T) 

roll 

( 8 ) 

roll sloshing 

<« s > 

rate gyro roll 

(0 r > 

swivel engine 

<*e> 


The control command equation, including transfer functions for all filters and 
the actuator characteristics, is the final equation. 

Two numerical methods are presented for solving the system's eigen- 
values: the characteristic equation and the matrix iteration approaches. In 
both cases, the system is solved for the normalized eigenvectors, which provide 
additional information about the system. 


GENERAL APPROACH 


Basic Equations 

The equations of motion, the control equation, and the equation describing 
the response characteristics of the rate gyro are homogeneous, linear differen- 
tial equations. These equations are transformed into a set of homogeneous, linear 
simultaneous algebraic equations by assuming solutions of time dependency in 
the form e st by which all differential quotients with respect to time are replaced 
by the complex operator. 


S = CT + ico . 

(1) 

With the coordinate (or unknowns) denoted as Xj(j = 1, 2, 3, . . . , n) and 
the coefficients of Xj as dy ^ g ) ( i, j = 1, 2, . . . , n) , the set of equations reads 

d lMs) X,+d 12(s) X * + - 

. . . . + d, . .X =0 
ln(s) n 

d 2 l( S > X ‘ +d 22( B ) X * + - 

. . . . +d 0 . .X =0 
2n( s) n 

(2) 

... ... 

.... ... — o 

d nl, s) X < + d „2(s) X * + - 

. . . . + d . .X = 0 
nn(s) n 
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To keep a consistent flow of data from other programs needed in equation 
(2) , the following order and definitions of the X. and equations are used 
throughout: ■* 


Number (i) 
1-2 

3 

4-5 

6 

7 

8 


Definition 
Torsional modes 
Roll 

Roll sloshing 
Rate gyro roll 
Swivel engine 
Control 


Symbol 


T. 

l 


£ si 


By using matrix notation, equation (2) can be more conveniently ex- 
pressed in the form 


D (s) { X.}=°, 

where all the elements of the matrix D. . except the nth row have the form 

(s) 

d... . = s 2 A.. + sB.. + C„ 

13 (s) ij ij ij 

(i = l — n - 1) (j = 1 — n) 


(3) 


(4) 


and the elements for the nth row, which describe the control equation and the 
filter characteristics, as 


d , =K.T. / . 
nj(s) j j(s) 

(j = 1 — n) . 


(5) 
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The transfer functions T.. , are of the form 

j(s) 


m n=0 

T = 

J 10 


E m" 

a=0 J 
10 

E y m 

v^— A J 


With the general equation (3) in this form, two basic approaches are 
available for determining the eignevalues [1]: (1) the matrix iteration and 
(2) the characteristic equation approach. (The characteristic equation is ob- 
tained by expanding the determinant and solving for its roots.) Once the roots 
( eigenvalues) are determined, it is important to solve for the eigenvector for 
additional information. This is done by assuming a value for the engine deflec- 
tion command angle (X or /3 ) and solving for the resulting eigenvector. 

n C 

Usually the eigenvectors are normalized to or /3 equal to one. 

Matrix Iteration for Obtaining Eigenvalue 
Equation (3) , written as 

D( S) {X.}=°, (7) 

states the transformed equation of the system in matrix form. The problem is 
to find nontrivial values of s and (X J that satisfy equation (3) . The eigenvalue 

problem is nonlinear because the coefficients of D . . are usually higher degree 

( S) 

polynomials in s. Starting with an approximate value of s and (X^) , it is possi- 
ble to find a set of linear nonhomogeneous equations for the correction terms 
that have to be added to the approximate value of s and (X.) . By developing 

equation ( 3) into a Taylor series at the point s and neglecting terms of higher 
order, the iteration procedure (1) is 


D. >(AX . + D' }{X }As 
( S n> J n_1 <s n > l n 


-D. ,X. =-Y. , 

(s n ) jn jn 
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where 




To obtain d 


(s) 

ds : 


each element of D, . is differentiated with respect 
(s) 


to s. The numerical procedure for equation (8) is as follows: 

Step I. - Insert first approximation of s = s ^ - obtained from approx- 
1 — n n 

imate root programs, sponsor load in choice, or the natural frequencies of the 

modes - into equation (3) and solve for the approximate eigenvectors (X.) by 

( 0 ) J 
setting equal to one. 

Step n . - Determine the right-hand side of equation ( 8) (-Y.) by using 
the eigenvectors obtained in Step I. ■* 


Step HI. 
be solved for. 

Step IV . 

Step V . 


- Solve equation ( 8) for As. 


The AX. are not needed, but could 
J 


- Calculate = s^ + A s . 

- Check s versus a constant (0.001 + i 0.001) . 


a. If either the real or imaginary part is larger than 0. 001, 

begin Step I again with ^ as s ^ and repeat the above 
step. 

b. If both real and imaginary parts are smaller than 0. 001, 
the program has converged to an eigenvalue. Using the 
eigenvalue, compute the ratios of the eigenvectors for the 
engine deflection (3^ as described previously. 

Step VI . - Print out results. 

Experience has shown that normalizing the eigenvectors in Step I by 
assuming X r equal to one does not always give the best convergence. The con- 
vergence problem can be relieved by selecting the X. for normalization to cor- 
respond with the eigenvalue that you wish to find. This means that provisions 
must be made in the program for interchanging rows and columns of equation ( 4) 
before starting the iteration procedure with Step I. 
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Characteristic Equation Approach 

The other approach for solving equation (4) for its eigenvalues is using 
the computer to expand the determinant of the coefficients into a poly- 

nomial. In general, two restrictions are imposed on solving the determinant: 

( 1) the determinant must not be larger than 10x10 (in our case, it is smaller) 
and (2) the elements must be single polynomials in s, not ratios of polynomials. 
Equation (4) does not usually meet both of the requirements and must be altered. 
First, the important modes must be determined (not necessary for our system) ; 
and second, the determinant must be expanded by minors along the last row, since 
the last row contains ratios of polynomials. 

Equation ( 4) may be rewritten as 


D 


(s) 


j=l 




- 0 , 


(9) 


where I?, . is the minor of the determinant D . 

(s) ( 

is a ratio of polynomials and can be wriiten as 


taken on the last row. 




N 


j(s) 

j(s) 


( 10 ) 


Equation ( 9) is rationalized and its numerator set equal to zero to solve 

the characteristic equation. The eigenvalues are the roots of the numerator. 

D . , is written in the form 
(s) 


D 


(s) 



(-l) n ' J D. j . N 

(s) j(s) 


K 

n p 

K*j 



0 . 


( 11 ) 


After the eigenvalues have been determined, equation ( 8) must be solved for 
the normalized eigenvectors with an option to bypass. 
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The matrix iteration has two advantages: (1) only important roots are 
obtained, and (2) the size of the system is not critically limited. The disad- 
vantages of this method are: (1) more machine time is required, and (2) con- 
vergence is poor for high gains with a gain root locus study. 


The characteristic equation approach is advantageous because all roots 
can be found, less machine time is required, and gain root locus studies offer 
no adverse convergence problems. However, this approach can solve only 
small systems. Systems with higher than 30th order create problems. 


Frequency Response 


the 

the 

the 

for 

set 


The frequency response of the system to a sinusoidal forcing function at 

control thrust vector point is easily obtained. From equation (3) we let <r of 

assumed solution s be equal to zero and assume various values for <x> within 

interested frequency region. For each value of w , equation (3) is solved 

the unknown X. by assuming X ( corresponding to /3 equal to one) . We 
J n c 

Kj equal to zero ( see definitions of elements) . 


Using these frequency-dependent unknowns, a structural transfer func- 
tion can be obtained between a force applied at the swivel point and any sensor 
location. 

1. Between the thrust vector and a position gyro 


F 


e (w) + y t m(^) e 

F £ F fi(x ) 
M g 


( 12 ) 


( x is the position gyro location) . 
g 

2. Between the thrust vector and a rate gyro, 


K(x ) 

r 

F 


= ico 




(13) 


Other frequency responses can be obtained from equation (13) by solving 
for the ratios of certain unknowns. 
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BASIC PROCEDURES FOR ANALYZING 


The analysis to be done is either a conventional gain root locus or a phase 
root locus with a constant gain. In conventional graphical techniques, the phase 
root locus is much more difficult to construct than the gain root locus; however, 
for the numerical procedures presented in the report, these problems are elim- 
inated. A combination phase and gain program provides more information than 
a conventional root locus diagram [2] . 

For a phase root locus equation ( 5) is written as 


K p 
n c 


n-1 

-A T KfT .X. 
(s) J J(s) J 


(14) 


Because this is the control equation written in general form with all the 

transfer functions of the various control loops, a phase root locus can now be 

computed for the total system or for each control loop separately by writing 

K?'s in the form 
J 


id. 

K? = K.e 3 
3 J 


and 


K 

n 



(15) 


Equation ( 13) then becomes 


e 



c 





j=i 


j(s) 


id. 


K.e 3 
J 


X. 

J 



(16) 


Now K contains the open loop gains of each loop as a parameter and is 
kept constant as the phase 6. is varied and the roots of equation (3) are com- 
puted for each value of 6. . 
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The matrix iteration procedure works much better for the phase root 
locus, since the complex elements introduced into the control equation do not 
cause any convergence problems. This is not always true in the root solving 
procedure for polynomials. 


EQUATIONS OF ELEMENTS 


The coefficients of s 2 A.. are as follows: 

U 

Torsional Modes ( i = 1 — 2) 

j -1-2* 

A.. =M .6.. 

13 ri 13 

j =3-6 

A.. =0 
U 

j =7 



= deSefl . . 

C(x t ) 
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j =8 


A i8 = ° 


Roll (i = 3) 

j =1 


A 33=° 


= 3 


A ~- + <^-W> + 2 M sn dS n 

n=l 


j = 4 — 5 
2 

* v at j n = 1 when j = 4 

A = > M ds „ , c 

3j LJ , sn n n = 2 when l = 5 
J n=l 

j =6 

A 36 = ® 


j =7 

A 37 = deSe 

j =8 

■A-38 = 0 

Sloshing Roll (i = 4 — 5) 
j =1-2 

A « = ° 

) =3 
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A = ds. 
i3 1 


j =4 — 5 


A.. = 0 when i # j 
ij J 


A.. = 1 when i = j 
U J 


i =6 — 8 


A.. =0 
ij 


Rate Gyro Roll (i =6) 
j =1-5 


A er° 


j =6 


Afie ~ U> 2 

r 


j =7 — 8 


A 6i'° 


Swivel Engine ( i = 7) 


j =1-2 

deSefl.. . 

J(x) 

A = — 5~ — 

7j 0 co 2 
e e 


j =3 


A _ deSe 

A?3 " e to 2 

e e 
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j =4-6 



3 =7 


_1_ 

A77 — 2 

CO 

e 

j = 8 
a 78 = 0 

The coefficients of s B.. are as follows: 

ij 

Torsional Modes (i = 1 — 2) 


j =1^2 


B.. = 2£ M 

1J T1 


.M .6.. 
ri ti ij 


j =3-8 

B.. = 0 
ij 

Roll (i = 3) 

j =1-8 

B 3j = ° 

Sloshing Roll ( i = 4 — 5) 

j =1-3 

B.. = 0 
i] 

j =4 — 5 

B.. = 2f oj . 6 .. 
ij si S1 L 3 
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j -6 — 8 


B.. = 0 
ij 


Rate Gyro Roll ( i = 6) 
j =1-2 


B^. =-0., 

6j j(x<£r) 


j = 3 

B G3 = “ 1 

j =4 — 5 


B er 0 


j =6 
2£ 


B GG ~ 


Ct) 


j =7-8 


B 6j = ° 


Swivel Engine (i = 7) 
j =1-6 


B,,. = 0 
7j 


j = 7 

2£ 


b 77 - 


Ct) 
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j = 8 

2£ e - 
B 78 = - ~~ 

00 

e 


The constant coefficients of C„ are as follows: 
Torsional Modes (i = 1 — 2 ) 


J 


i — 2 


C . = co ?M . 6. . 

1] T1 T1 1] 

j = 3 — ► 6 


C.. =0 
1J 


j 

C.„ = 8 , v deF 
i 7 M(x t ) s 

j =8 


C i8 


Roll (i = 3 ) 

j =1 


° 3 j 


j =7 

C37 = deF 


j =8 
C38 = 0 


14 


Sloshing Roll ( i = 4 — 5) 


j =1 -3 


C.. =0 
ij 


j =4 — 5 


C.. 

i] 


= w 


. 2 5 .. 

si 1J 


j =6-8 

C.. = 0 
ij 

Rate Gyro Roll ( i = 6) 


j =1-5 


C 6J = ° 


j =6 

C 6G = 1 
j =7 — 8 


C 6i=° 


Swivel Engine ( i = 7) 


j =1-6 


C 


7j 


= 0 


3 = 7 

C77 = 1 

3 =8 

C78 = "^1 
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Control Equation (i = 8) 


j =1-2 


1U * 

C_. = — a e. A A T K<e 1 
8j or j(x<£r) (s) g(s) 1 


C33 = -a AT. I^e 
or (s) g(s) 


<v° 


c 86 --a lr A (g) T r(s) K 2 e 


j =7 


C07 — 0 


j =8 


Coo — 
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APPENDIX 


Roll Equations 


The equations of motion for the roll stability program are derived using 
Lagrange's equation. The torsion, roll, roll sloshing, and swivel engine equations 
were derived, while the rate gyro equation was obtained from Reference 3. The 
longitudinal acceleration was eliminated from the equations of motion because of 
its negligible influence on the final results. 

1. Lagrange's Equation 



where L = T - V. 

T is the kinetic energy, and V is the potential energy of the system, while 
D is the dissipation function, and Qj's are the generalized forces. 

The generalized coordinates, q., for the derivation are: 



Torsional twist 

q 2 = o 

Rigid body roll angle 

q s = £ s 

Roll sloshing displacement 

q4 = 0 r 

Rate gyro roll 

q 5 = 0 e 

Swivel engine 

2. Kinetic Energy 

a. 

Airframe 


The kinetic energy of the airframe depends on the rotation of each 
element, which is the sum of the rigid body roll angle and the torsional twist. The 
angular velocity about the missile centerline is 
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(2) 


A(x) 


= ^Z 


T 9 , . . 

M M(xJ 


The kinetic energy is given by 


T 


A 




(3) 


where I is the axial distribution of the roll mass moment of inertia, 
(x) 


= i fl 

2 i <*> 


0 2 + 


29 ^ T K 8 M<x,)j 


dx 


+ 




dx 


T = - 
A 2 






(x) 0 M(xt) dX 



(x) 




dx . 


But / I (x) dX = I R ’ < 4) 

A 

where I R is the total moment of inertia of the airframe about the centerline and 

/ I (x)%(x t ) dX= ° ( end torc l ue ) • ( 5 ) 
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Also, 


J I (x) 0 i(x t ) 0 j(x t ) dx 


1J \ i = i = M Ti 


( 6 ) 


where 5 . . is the kroneeker delta. 

13 

Therefore, 


2 . ± 


T A = 2 V + 


2jM t 2 


(7) 


b. Liquid propellant 

The outer clustered tanks of a Saturn type space vehicle are com- 
bined into two tanks for the spring-mass model (Fig. 1). Springs and dashpots 
are used to attach the oscillating masses to the tank wall. 

The angular velocity at any station x is 


V 


s(x) 


= ds 0 
n 


+ k 


( 8 ) 


where ds^ is the distance from the centerline of the vehicle to the sloshing mass 
(n = 1,2). 

Hence, the kinetic energy is given by 


T v = ~ M V . 2 = — Yj M < ds 0 + 1 ) 2 dx 
(s) 2 u sn s(x) 2 u sn n ^s 

n n 


t = 4 Y 'm ds 2 0 2 + Vm ds e£ 

(s) 2 ^ sn n ^ sn n S s 2 u sn^s 

n n 


(9) 
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c. Swivel engine 


The angular velocity of the gimbal rotation is given by (Fig. 2): 


V 0e(x) ° + ) 


Therefore, the linear velocity at the gimbal station is 


V = - de ( 0* + YjT 6 . . ) , 

e(x) l ^ M M( x t ) / 


where de is the distance from the missile centerline to the swivel engine. The 
velocity of a point on the engine relative to the gimbal is given by We' Con- 
sequently, the velocity of a point on the engine is given by 


V . - d. 8 + de I 0 + Yi 6 . . , 

e(x) (xfe 1 ^ n n (x t ) j 


and the kinetic energy of the system is 


r =4 f M , v d. J + de 6 + YjT 6 . J dx, 
e 2 J e(x) (x) p e l ~ ix M( x t )J 


But f M . d 2 . dx = 0 
J e(x) (x) 


when 0 is the moment of inertia of the engine about the gimbal plane; and 

G 


de! / M e(x) dx =I r=°- 


where I is the moment of inertia of the engine about the missile centerline. 
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f M d dx = Se , 
J e(x) (x) 

tJ 


( 16 ) 


where Se is the static moment of inertia of the engine about the gimbal. 
Hence, the kinetic energy of the engine is 


T = f 6 ji 2 + deSe0 p + deSe/3 + Y. t 6 . , . 
e 2 e'e K e u n (x ) 

M * 


(17) 


3. Potential Energy 
a. Airframe 


airframe. 


The torsional vibration comprises the total potential energy of the 


V = - Im u> 2 t 2 , 
A 2 u t n th n 

/Li 


(18) 


where w 
mass. 


is the torsional natural frequency and M 

T/i 


is the generalized torsional 


b. Liquid propellant 

The potential energy usually arises from three sources, but in roll 
only the strain energy contained in the spring of the mechanical analogy is used. 


V = - k 
s 2 n*s 




sn sn*s 


(19) 


n 

where k is the spring constant of the slosh mass, 
n 

c. Swivel engine 

The swivel engine is coupled with the roll control system, and the 
dynamics of the swivel engine can be described by the simple spring-mass-damper 
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system. The actual potential energy is determined by an effective spring constant 
and damping coefficient of the gimbaling system. 


V e-K<V' , c )l= 2 


9 or (/3 
e e e 


/U ! 


( 20 ) 


where k is the spring constant and w is the natural frequency of the swivel engine, 
e e 

4. Dissipation Function 

a. Airframe 

The dissipation of the empty airframe arises from its structural 
damping, which is a function of the torsional twist and not of its frequency. 


D = ~ w g T 2 , 

A 2 u TjU Tp Tp p 


( 21 ) 


where g is the torsional damping of the structure. 

Tp 

b. Liquid propellant 


The dissipation of the liquid propellant is represented by an equiva- 
lent dashpot associated with the mechanical model. 


D 

s 


1 ? 2 

= — C £ 

2 s^s 


=ii> 


n 


w g , 
sn sn sn s 


(22) 


where g gn is the propellant damping, 

c. Swivel engine 

The damping of the swivel engine is of minor importance and can be 
neglected, but the viscous damping of the actuator is proportional to the difference 
in the velocities of the command signal and the engine deflection; 


D, 




( 23 ) 


where g g is the swivel engine damping. 
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5. External Forces 


The only external force considered in this derivation is the swivel thrust. 
The drag, aerodynamic forces, and forces caused by the flowing propellants are 
neglected. 


Q = deF 0 . sin ( -R ) = - deF 0 , , sin B , 

t/u s M(* t ) e s /x(x t ) p e* 


and we may write 


sin R 

e 




= -deF 0 J3 
s /i(x t ) e 


(24) 


Q e = deF s sin(-3 e ) = -deF^ 


(25) 




= 0 . 


(26) 


6. Equations of Motion 

a. The kinetic potential, L, is given by 


L=-I0 2 + -£m t 2 + - ds 2 0 2 +Y.M ds 0 | 

2 R 2 *-• tu u 2 u s a n u sn n S s 
H n n 


+ 77 Xm I 2 +T0 /? + deSe0/3 + deSe/3 T.r 0 , . 
2 ^ sn b s 2 e e 'e K e u /j. ^(x^.) 


“ w2 t 2 "I" Xm a) 2 | 2 0 u) 2 ^ 2 

2 ^ T/i tu ju 2 “ sn sn*s 2 e e*e 


n 


+ 0^UJ 2 /? /3 - 4 0 w 2 ^ 2 . 
e & e c 2 e e K c 


( 27 ) 
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b. The dissipation function, D, is given by 


D = -Y,M w g t 2 +-tY,M io g | 
2 ^ tju T/u T/u m 2 L-t sn sn sn s 
u ^ n 


+ -g0wj3 2 -g0ci>/5/3 +-^g0ct)/3 2 
2 s e e e p e 6 e e e'Vc 2 °e e e K c 


( 28 ) 


c. The equations of motion are written with the aid of Lagrange's 

equation: 


_d 

dt 


\ 

9L _9L _9D 

8q. ’ 9q. 9q. ^i * 

l/i l 


( 1) Torsion equation ( t ) 


= Ym t + deSe/5 
8t ^ tj ix /x e 




2 * 

M 


M(x,) 


_d 

dt 


9L 

ar 


= Y m t + deSe0 

u TH n m(x ) r e 

p 1 



Ym cu 2 t 

T/X T/Ll jU 


8D 

af 

ju 


Z M 

/u 


W g T 
T/X T JU T/LX /X 


Q = - deF 0 . . fl 

r/x s /x(x,) r e 
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Hence, the torsional equation is 


Z M t + TIm u> g t + y,M w 2 t + deSe0 

Til II T/l T/l T/f II 1 Til Til II 


TH n 


TH Til T/X IX 




TfX TfX fX 


V( x t )' S e 


+ deF 0 . .3 =0 

s n(x re 


( 29 ) 


(2) Roll equation (0) 


9L 

90 


= L0 +Xm ds 2 0 + V M ds £ + deSeS 
K " sn n u sn s 
n n 


d 9L 1 •• v o •• v' *• •• 

tt rr =L,0 +X M ds 2 0 +Xm ds | + deSe/J 
dt l 90 R ^ sn n ^ sn n*s *e 


n 


n 


9L 

90 


= 0 


9D 

90 * 


= 0 


Q. = - de F 0 
0 s e 


Hence, the roll equation is 


L.0 + y.M ds 2 0 +y,M ds £ +deF B + deF 5=0 
R u sn n u sn n*s s'e s p e 

n n 


(30) 
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(3) Sloshing roll equation U g ) 


-^ = T> ds 6 + 5> i t 

Lj sn n ^ sn^s 


’s n 


n 


_d_ 

dt 


v 

\ d(s l 


=t.m ds 0 + § 

Lj sn n LJ sn s 
n n 


-£m u> j | 

9£ ^ sn sirs 

n 


=y M u g £ 

d£ u sn sn sn’s 
S s n 


Q {s =0 


Hence, the sloshing roll equation is 


V m ds fl + T M £ + Tm uj g £ + T> | — 0 
ZiVn ^ sn 5 s ^ sn sn B sn 5 s ^ sn sn s s 

n n n n 


Dividing through by the sloshing mass, 


ds 0 + £ + 2w t £ + co 2 £ = 0. 
n ? s sn b sn s s sm“s 


(.31) 


(4) Rate gyro roll equation (0 r ) 

The rate gyro equation is obtained from Reference 3, but can 
easily be converted to roll. 


26 


The rate gyro is included as a separate equation because it is nonideal. 
Hence, the rate gyro roll equation is 


9 2£ 

r i 

CO 2 + CO 


9 + 9 = 
r r 


9 + 


£t0 

LJ a 


H n(x<pr)j 


( 32 ) 


9 

- Y t 9 - 9 + — i - 

M M(x^r) w 2 r 



9+9=0. 
r r 


(32 a) 


(5) Swivel engine equation (p ) 

c 


3L • ♦ V' • 

— r- = 9 B + deSe0 +deSe ),t 9 . . 
dp ere u n /i(x.) 

K e n t 


_d_ 

dt 



= 9 p 
e^ e 


+ deSe0 + 


deSe T. t 0 . v 
" M M(xJ 


= 9 u 2 p - 9 co 2 /3 
9/3 e e H c e e p e 



= g 9 B 
& e e e p e 


- g 9 on p 
e e ec 



= 0 


Hence, the swivel engine equation is 
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deSe^ 


• • • * 


r 0 . + deSefl + 0 / 3 + 2£ u> 0/3 + 0 w“/3 

u ulxj e'e e e e'e e e e 

H t 


- 2£ u> 0 B ~ 0 w 2 /3 = 0 . 
e e e c e e c 


Dividing by 0 or : 

6 G 


deSe Y)t 0 , , *• 0 . 

m ^ ^ X t deSe g* ^e 


2f, 


0 W* 

e e 


+ — 5” + /3 + B - /3 - /3 — 0 . ( 33) 


? or ur w e e u> * c c 
e e e e e 


(6) Control Equation 


a or 0 + ^ T /n 0 ju(x0r)J + a ir 0 r 


(34) 


-a 7 !t 0 . , . - a 0 - a 0 + B =0 
or “ fx n(x<t>r) or lr r c 


(35) 


With filters and actuator 


a or ^ T |it%(x0r) A (s) T g(s) K le a or T g(s) A (s) 0K i e 


i6 e i6 A 

-a T ,A 0 K,e 2 + — — fi = 0 . 

lr r (s) (s) r 2 K c 

A 


(35) 
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FIGURE 1. SLOSH MECHANICAL MODEL 
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FIGURE 2. SWIVEL ENGINE MODEL 
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